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INTRODUCTION 

It is important to accurately predict the pressure, 
temperature, as well as the amount of 0-ring erosion in 
the space shuttle Reusable Solid Rocket Motor 
(RSRM) joints in the event of a leak path. The 
scenarios considered are typically hot combustion gas 
rapid pressurization events of small volumes through 
narrow, restricted flow paths. The ideal method for this 
prediction is a transient three-dimensional 

computational fluid dynamics calculation with 

computational domain including both the combustion 
gas and the surrounding solid regions. However, this 
has not yet been demonstrated to be economical for this 
application due to the enormous amount of computer 
time and memory required. Consequently, all CFD 
applications in RSRM joints*' are steady-state 
simulations with solid regions being excluded from the 
computational domain by either assuming a constant 
wall temperature or assuming no heat transfer between 
the hot combustion gas and cool solid walls. 

Currently there are two computer codes, known to the 
authors, available to model the gas dynamics, heat 
transfer, and 0-ring erosion in the RSRM joint 
pressurization process. One is ORING2'^ which was 
developed at Thiokol Propulsion, and the other is JPR . 
which was developed at NASA Marshall Space Flight 
Center. A way to improve the current prediction 
technique is to modify the transient compressible flow 
calculation since the pressure, temperature, and 
velocity of the combustion gas are not calculated from 
the time-dependent Navier-Stokes equations. Instead, 
some empirical correlations are used to predict the gas 
temperature, mass flow rate, and other flow properties 
by assuming a quasi-steady state flow in a constant 
cross-section area pipe (i.e., there is no grid in the 
paths). Furthermore, ORING2 can only handle 
configurations with two volumes and two paths while it 
takes significant coding for JPR to do complicated 
configurations with more than two volumes. 


A new thermal-flow simulation code, called SFLOW, 
has been developed to model the gas dynamics, heat 
transfer, as well as 0-ring and tlow path erosion inside 
the space shuttle RSRM joints. The details are 
discussed in this paper. The SFLOW methodology 
eliminates some of the approximations inherent in other 
simulation and prediction tools. This is accomplished 
by combining SINDA/G® (Network Analysis, Inc.*), a 
commercial thermal analyzer, and SHARP®^'*^, a 
general-purpose CFD code developed at Thiokol 
Propulsion. The pressure, temperature, and velocity of 
the combustion gas in the leak paths are calculated in 
SHARP' by solving the time-dependent Navier-Stokes 
equations while the heat conduction in the solid is 
modeled by SINDA/G®. The two codes are coupled by 
the heat transfer at the solid-gas interface. The number 
of flow paths and volumes in SFLOW is limited only 
by the memory of the computer used to run SFLOW. 

Although SHARP® can solve one-dimensional, two- 
dimensional, as well as three-dimensional flow 
problems, the tlow inside paths is assumed to be one- 
dimensional in the current version of SFLOW to reduce 
the CPU time and memory requirements. This is a 
reasonable approximation due to the fact that leak paths 
arc usually narrow. The solid in SFLOW, however, can 
be one-dimensional, two-dimensional, as well as three- 
dimensional since it is not always a good 
approximation to assume the heat conduction in the 
solid region is one-dimensional, especially when the 
wall material is metal. This is feasible in terms of CPU 
time and memory because there is only one equation 
(i.e., conservation of energy) to solve in the solid 
region compared to five equations (i.e., conservation of 
mass, momentum, and energy) needed in the gas region 
if the flow is modeled as three-dimensional. 
Furthermore, a larger time step can be used in the solid 
calculation than that used in the gas calculation since 
the solid temperature usually changes much slower than 
the flow properties. 
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The SHARP® main program is converted into a 
subroutine so that it can be called from SINDA/G®. 
The input for this subroutine includes the heat flux 
from gas to wall, friction factor of the flow path, mass 
addition due to erosion, gas properties, grid, as well as 
boundary and initial conditions while the output is the 
pressure, temperature, and velocity of the gas for each 
flow cell at a specific time step. In SFLOW, the flow 
calculation (i.e., SHARP®) and the solid calculation 
(i.e., SINDA/G®) are decoupled from each other such 
that a smaller time step can be applied in SHARP® than 
that used in SINDA/G®. Furthermore, SHARP® can use 
more ceils for the flow solution than the number passed 
from SINDA/G® solid surfaces. 

As a general-purpose CFD code, SHARP® does not 
have the friction term in the governing equations since 
it is typically used for two-dimensional and three- 
dimensional flow simulations where the friction is 
implicitly taken into account by the viscous torce in the 
resolved near-wall region or by a wall function. In one- 
diniensional flows, however, the friction term has to be 
explicitly added due to the fact that the velocity 
gradient in the wall-normal direction does not exist. 
Similarly, a heat transfer term is added to the SHARP® 
equations because the thermal boundary layer is not 
simulated in SFLOW. Moreover, a mass addition term 
is added in SHARP® since the erosion of the wall 
material will generate mass. Finally, minor loss terms 
such as those due to sudden expansion or contraction 
and flow direction change are accounted for by 
specifying a loss coefficient in the SFLOW input file at 
the appropriate flow cells. 

The details of the solution scheme, including the 
modeling of gas dynamics, heat transfer, as well as O- 
riniz and path erosion, are discussed in the next section, 
followed by comparison of SFLOW predictions to 
exact solutions or experimental data. The test cases 
included Fanno flow where friction is important, 
Rayleigh flow where heat transfer between gas and 
solid is important, flow with mass addition due to the 
erosion of the solid wall, transient volume venting 
process, as well as some transient one-dimensional 
flows with analytical solutions derived by Cai . In 
addition, SFLOW has been applied to model the RSRM 
nozzle joint 4 subscale hot-flow tests'^, which simulate 
flows to the primary and secondary 0-rings. The 
predicted pressure, temperature (both gas and solid), 
and 0-ring erosion from SFLOW are compared with 
the measured data in this paper. 


GAS DYNAMICS AND THERMAL 
MODELING 


GAS DYNAMICS MODELING 

In SFLOW, the gas can be either in a flow path or a 
volume (i.e., cavity), which are treated very differently. 
The gas in a volume is assumed to be in quasi- 
equilibrium with uniform pressure, temperature, and no 
velocity whereas that in a path is solved from the first 
principles (i.e., conservation of mass, momentum and 
energy). 


Gas flow in paths 


The transient compressible flow in a path is modeled 
using SHARP®, which is a general-purpose CFD code. 
SHARP® solves the Reynolds averaged Navier-Stokes 
equations in one-dimensional form as 

dQ ^ 3(£,-£. ) _^ (1) 

dt dx 


where the unknowns are 
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( 2 ) 


In equation (2), A is the cross-section area, which is a 
function of both space .r and time /, and p and u are 
the Reynolds averaged density and velocity, 
respectively. The total energy is 

e = pj^c,r + Lrj (3) 


where T is the Reynolds averaged temperature and 
is the specific heat at constant volume. The inviscid 
flux term is given by 
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while the viscous term is 
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The total stress and heat flux include both laminar parts 
and turbulence parts as 




(6) 
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where ^ '■ and Pr‘ are the laminar viscosity and Prandtl 

number while Pr’^ is the turbulent Prandtl number. The 
turbulent viscosity ^ is zero for laminar flow while it 

is obtained by the widely used k—£ model for 
turbulent flow. The source term in equation (1) is 
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where th is the mass addition rate due to path erosion,/ 
is the Darcy friction factor, D/, is the hydraulic diameter 
of the path, K is the minor loss coefficient and q is the 

heat transfer rate per unit volume from the gas to the 
solid wall. The boundary conditions for the path 
flowfleld are obtained from the pressure and 
temperature of the volumes. A path has to be connected 
with a volume at one end; the other end can be 
connected to a volume or a solid wall, which could be 
either adiabatic or conducting heal away to the solid 
region. 

Note that, similar to other general-purpose CFD codes, 
SHARP® does not have the friction term shown in 
equation (8) since it is typically used for 2D and 3D 
flow simulations where the friction is implicitly taken 
into account by the viscous force in the resolved near- 
wall region or by the wall function applied. In one- 
dimensional flows, however, the friction term has to be 
explicitly accounted for due to the fact that the velocity 
gradient 3w/3ydoes not exist. Similarly, the heat 

transfer term in equation (8) is added to SHARP® due 
to the fact that the thermal boundary layer is not 
simulated in SFLOW. The mass addition terms are also 
added since the erosion or decomposition of the walls 
will generate this effect. Finally, minor loss terms such 
as those due to sudden expansion or contraction and 
turns or bends in the flow path are added in equation 
( 8 ). 


Friction factor in the flow paths 

The friction factor in equation (8) is obtained from the 
empirical correlation of Idelchik*^, which depends on 
both the shape of the path (i.e., circular, rectangular, or 
triangular) and whether the flow is laminar, turbulent, 
or in transition. 

According to Idelchik*^, the gas flow is divided into 
laminar, turbulent, and transitional regimes depending 
on the cutoff Reynolds numbers defined as 

Re, = 754exp(().0065D, ! e) (9) 
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where £ is the roughness of the flow path. The flow is 
laminar if the Reynolds number is below Rea, turbulent 
if the Reynolds number is above RCc , and in the 
transitional regime if the Reynolds number is between 
RCa and RCc. If there are two transitional zones, Rcb is 
used to determine which of the two zones the flow is in. 


The friction factor in the flow path is then determined 
based on the flow Reynolds number as: 

• for Re < Re, 



• for Re , < Re < Re, 


/ =4.4Re'^‘'"'exp| 
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• for Rc/, < Re < Re^ 
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for Rc„ < Re < Re< and Re^ < Re„ there is only 
one transitional zone 
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/ = (7.244 Re ■'’'^’-0.32)exp (17) 

(- 0.00 1 7- (Re^ - Re)- ) + 0.032 

Pressure and temperature in volumes 

Once the flowfield in the path is solved by SHARP®, 
the pressure and temperature in the volumes can be 
obtained from mass conservation 

(18) 

dt 

and energy conservation 

^{mcj)=ymh-Q (19) 

dt ^ 

where the summation is for all paths which connect to 
this volume, rii and h are the mass flow rate and 
enthalpy at the end of the path, is the rate of mass 
addition to the gas due to surface erosion, Q is the heat 

transfer rate from the gas to the solid boundary which 
includes the convective heat transfer as well as the heal 
transfer due to erosion, m, P, and T are the mass, 
pressure and temperature of the gas in this volume, 
respectively. In addition to equations (18) and (19), the 
ideal gas law 

pV = mRT (20) 

where V is the volume of the cavity, was used to solve 
the pressure p, temperature 7, and mass m of the 
volume. 

HEAT TRANSFER MODELING 

The convective heat transfer between the gas and the 
solid wall is modeled as 

q = hA,(T^,-TJ (21) 

where h is the heat transfer coefficient, is the 

surface area, and are the temperature of the gas 

and solid wall, respectively. This heat transfer rate is 
used in both SHARP® and SINDA/G® so that the total 
energy in the system is conserved. 

Heat transfer in paths 

The heat transfer coefficient in flow paths can be 
obtained from the Nusselt number as 

h = N^- ( 22 ) 

where k is the thermal conductivity. The Nusselt 
number depends on both the cross-section shape of the 


path and the flow regime. If the flow is laminar and the 
path is circular 

N„=436 (23) 

while for rectangular paths 

= i , 1 8 1 35 + 2.30595r”"''^^"^ (24) 


where 
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with Cl and h being the width and height. For turbulent 
flow, the Nusselt number is calculated using the 
following empirical correlation 


=- 


/ Pr Rc 


^(Pr'''.(l.07-l2.7V//8{P^-''-lj 




(26) 


where / is the friction factor calculated from equation 
(12) through (17), and fig and are the viscosity 

evaluated at the average gas temperature and wall 
temperature, respectively. In the transitional regime, a 
linear interpolation between the laminar and turbulent 
Nusselt number is applied. 


Jet impingement heat transfer 

The jet impingement heat transfer correlation used in 
SFLOW is the same as that in ORING2 and JPR, which 
depends on the standoff distance to diameter ratio as 
well as whether the flow is laminar, turbulent, or in the 
transitional regime. The heat transfer coefficient is 
obtained from the Stanton number as 


h = S,c — 
' '■ A 


If the flow is laminar, the Stanton number is 

xl/6 


S. = 0.763 — 

'Re 


1 (TA 


(27) 


(28) 


where the Reynolds number Re is calculated at the jet 
exit, is the gas temperature at the jet exit, and 7^, is 
the temperature of the wall surface. For turbulent flow 
with a standoff distance to jet diameter ratio L! D <2.6 
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while for Li D >2.6 
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5, = 0.442 Pr'“" 1 + 
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(30) 


where Le is the Lewis number. The velocity gradient in 
the above equations is 

V',™. =1 (31) 


for L! D< 3.4, 




(32) 


for 3.4<L/D<8.4,and 
1 -expl 


V = - 

^rad 


0.2315L/D-0.74 


0.13L/D-0.39 


(33) 


for Lf D> 8.4. If the flow is in transitional regime, the 
heat transfer coefficient is obtained by linear 
interpolation between laminar and turbulent regimes. 

Heat transfer in volumes 

The heat transfer from the gas in a volume to the solid 
boundary can be modeled in four different ways 

• Using the impingement jet heat transfer correlation 
described above 

• Using the heat transfer coefficient in the paths 
connected to this volume 

• Using a conduction length as h-k/l 

• Using a user-specified heat transfer coefficient 

The user of SFLOW specifies which of these methods 
should be applied to calculate heat transfer coefficient 
for all the gas-solid interfaces in all volumes. 

EROSION MODELING 

The erosion model used in SFLOW is the same as that 
in ORING2 and JPR. Specifically, the erosion rate is a 
function of heat transfer coefficient between gas and 
solid, the gas temperature, as well as the wall 
temperature. Erosion has the following effects on the 
gas: 

• Increases the cross-section area of the path or the 
volume of the cavity 

• Adds mass to the gas 

• Adds energy to the gas 


VALIDATION CASES 

For all the validation cases shown in this section, there 
are only two volumes connected by one flow path. The 
pressure and temperature in one or both volumes are 
specified as input whereas the pressure, temperature, 
and velocity in the flow path are calculated using 
SFLOW. Most of these tests are for SHARP® in 
solving one-dimensional How problems with friction, 
heat transfer, mass addition, and area change since no 
modification is made to the commercial thermal code 
SINDA/G 

FANNQ FLOW 

For air entering an adiabatic lOO-ft-diameter duct with 
a total pressure of 3183.65 Ibf/fr and total temperature 
of 540°R, it can be shown analytically that the inlet and 
outlet Mach number will be 0.8 and 0.9, respectively, if 
the friction coefficient is assumed to be 0.0578, the 
pipe length is 100 ft, and the outlet pressure is 1829.12 
Ibf/fr. SFLOW was used to simulate this test case and 
the results are compared with the analytical solutions in 
Table 1. It is clear that the error in the predicted Mach 
number is smaller when more How cells are applied. 
However, the error drops much more from 20 cells to 
200 cells than that from 200 cells to 1,000 cells. For 
this particular case, 200 cells are enough to keep the 
error in both inlet and outlet Mach number below 
0.55%. 


Table 1, The Mach Number Predicted by 
SFLOW for the Fanno Flow Test Case 



20 cells 

200 cells 

1,000 cells 


0.8232 

0.8043 

0.8027 

error in 

2,90% 

0.54% 

0.34% 


0.8900 

0.8954 

0.8958 

error in 

1.11% 

0.51% 

0.47% 


RAYLEIGH FLOW 

For air entering a frictionless 0. l-ft-diameter, l(X)-ft 
long duct with a total pressure of 1481 Ibf/fr and total 
temperature of 524°R, it can be shown analytically that 
the inlet and outlet Mach number will be 0.2 and 0.25, 
respectively, if the heat addition is assumed to be 
60.367BTU/lbm, and the outlet pressure is 1389 Ibf/ft^. 
SFLOW was used to simulate this test case and the 
results are shown in Table 2. The error in the predicted 
inlet Mach number drops more from 20 cells to 200 
cells than from 200 cells to 1,000 cells. For the outlet 
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Mach number, the error actually increases from 20 cells 
to 200 cells, but decreases slightly from 200 cells to 
1,000 cells. For this particular case, 200 cells are 
enough to keep the error in both inlet and outlet Mach 
number below 0.65%. 


Table 2. The Mach Number Predicted by 
SFLOW for the Rayleigh Flow Test Case 



20 cells 

200 cells 

1 ,000 ceils 


0.1954 

0.1987 

0.1991 

error in Af 

2.30% 

0.65% 

0.45% 

A/„„, 

0.2498 

0.2484 

0.2485 

error in ^ out 

0.08% 

0.64% 

0.60% 


p = const 
p = const 

2C„(C„-t + C) (35) 

2C,;r + C, 

if both friction and heat transfer are neglected. This 
transient How with area change case is simulated using 
SFLOW with a uniform grid of 20 flow cells by 
assuming C^^ = C^ = = \ ‘^nd = 100 . The flow 

path is from x = 0 to X = 20 m and the time is from 
r = Osec to r = I sec. The velocity at the inlet and outlet 
from SFLOW are compared with the exact solution in 
Figure I, which shows a very good agreement even 
though only 20 flow cells are used. 


MASS ADDITION 

For air entering a frictionless adiabatic 0. 1-ft-diameter 
duct with a total pressure of 1708 Ibf/ft' and total 
temperature of 525°R, it can be shown analytically that 
the inlet and outlet Mach number will be 0.5 and 0.6, 
respectively, if 0.021 Ibm/sec of air at a temperature of 
525"^R is added to the flow and the outlet pressure is 
142 1 Ibf/ftl SFLOW was used to simulate this test case 
and the results are shown in Table 3. With 200 or more 
flow cells, the enor in the inlet Mach number is less 
than L4% while that in the outlet Mach number is less 
than 0.35%. 


Table 3. The Mach Number Predicted by 
SFLOW for the Mass Addition Test Case 



20 cells 

200 cells 

1,000 cells 

A/,„ 

0.5088 

0.5068 

0.5069 

error in 

1.76% 

1.36% 

1.38% 


0.5917 

0.5979 

0.5984 

error in M „ut 

1.38% 

0.35% 

0.27% 


TRANSIENT FT.OW WITH AREA CHANGE 

All the test cases shown above are steady state 
problems so only the SFLOW predictions at long times 
were compared with the exact steady state solution. In 
this and the following three sections, SFLOW was 
tested using one-dimensional unsteady flow cases with 
analytical solutions derived by Cai*^. 

It can be shown that for one-dimensional compressible 
flow in a circular pipe with a cross-section area 






(34) 


the exact solution of the governing equation is 


Transient Flow With Area Change 



Figure 1 . Comparison of the Predicted 
Velocity From SFLOW With the 
Analytical Solution of Equation (35) 


TRANSIENT H.OW WITH HEAT TRANSFER 

It can be shown that for one-dimensional compressible 
flow in a circular pipe with a heat transfer rate per unit 
mass of 


Y 

C,.t + C, 

y-1 C, 1 

[c,f + c, 


(36) 


the exact solution of the governing equation is 
p = const 




P = 


C.(Q + C,)' 

(C,.t + 

C..C + C, 

U = 


(37) 


Cj + C\ 
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if the cross-section area is constant and the friction is 
neglected. This transient flow with heat transfer case 
was simulated using SFLOW with a uniform grid of 
200 How cells by assuming C, = C 4 = 1 • 

= C 5 =5000 and C = 15. The flow path is from 

jc = 0 to JC = 20 m and the time is from f = 0 sec to 
t = 1 sec. The velocity and temperature at the inlet and 
outlet from SFLOW are compared with the exact 
solution in Figure 2 and Figure 3. The velocity profile 
at both the inlet and outlet agree very well with the 
analytical solution and the agreement for the 
temperature is also reasonable. Although not shown 
here, the results using 20 cells are much worse than 
those shown in Figure 2 and Figure 3 and it is expected 


Transient Flow With Heat Transfer 



Figure 2. Comparison of the Predicted 
Velocity From SFLOW With the 
Analytical Solution of Equation (37) 


that a better temperature prediction would be obtained 
by using even more flow cells. 

TRANSIENT FLOW WITH FRICTION AND HEAT 
TRANSFER 

It can be shown that for one-dimensional compressible 
flow in a circular pipe with a heat transfer rate per unit 
mass of fluid 




16A 

;r/V 


(38) 


the exact solution of the governing equation is 

p — const 
p = const 

aJa 

u= ^ - 


(39) 


if both the cross-section area and friction factor are 
constant. This transient flow with heat transfer and 
constant friction factor case was simulated using 
SFLOW with a uniform grid of 20 flow ceils 
and /= 0.004/ The flow path is from ^ = 0 to 
jc =20 m and the time is from / = 5 sec to r = l 0 sec. 
The velocity at the inlet and outlet, which are the same 
at any given time according to equation (39), , from 
SFLOW are compared with the exact solution in Figure 
4 , which shows a very good agreement even though 
only 20 flow cells are used. 

Transient Flow With Friction 
and Heat Transfer 


Transient Flow With Heat Transfer 



Time (sec) 

Figure 3. Comparison of the Predicted 
Temperature From SFLOW With the 
Analytical Solution of Equation (37) 



Time (sec) 

Figure 4. Comparison of the Predicted 
Velocity From SFLOW With the 
Analytical Solution of Equation (39) 
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Velocity (m/s) 


TRANSIENT FLOW WITH AREA CHANGH. 
FRICTION AND HEAT TRANSFER 

II can be shown that for one-dimensional compressible 
flow in a circular pipe with area/4, friction factor/, and 
heat transfer rate q given by 


exp(-^;r/C,jc)-i- in 

4 

(40) 

(Cj exp(- / C, .v)+ /C, In) 

(41) 

^ ( C, (c, exp(- V^/C,.t)+ 

(42) 

Vc,[ C,; + C, 


the exact solution of the governing equation is 

p = const 
p = const 

_ C;(c,exp(-77r/C|jr)+-y/C^/7r) ( 43 ) 

C^t + 


This transient flow with area change, friction, and heat 
transfer case is simulated using SFLOW with a uniform 
grid of 20 flow cells by assuming C^^ =C^ =C^ = U 

C| = 1000 and c, =15- The flow path is from a: = 0 
to X = 20 m and the time is from t = 0 sec to 
/ = 0. 1 sec. The velocity at the inlet and outlet from 
SFLOW are compared with the exact solution in Figure 
5, which shows a very good agreement even though 
only 20 flow cells are used. 

Transient Flow With Area Change, 
Friction and Heat Transfer 



Time (sec) 

Figure 5. Comparison of the Predicted 
Velocity From SFLOW With the 
Analytical Solution of Equation (43) 


VOLUME VENTING 

The cases discussed above focus on the gas flow in 
paths since only pressure, temperature, and velocity of 
the gas in the pipe were compared with the analytical 
solutions. To validate the volume pressure and 
temperature algorithm, the SFLOW prediction in a 
volume-venting experiment is compared with the 
measured data. In the experiment, a tank with a volume 
of 62.024 in^ is connected to a 2.27 in circular pipe 
with a diameter of 0.072 in, which is open to ambient 
conditions at 12.5 psia. The tank is at ambient 
condition initially and the tank pressure is increased to 
about 258 psia from f = — lOsec to / = 0 with the 
valve between the lank and pipe closed. This valve is 
opened at / = 0 and the tank pressure begins to fall. 
The SFLOW prediction starts from f = 0 using 20 
flow cells in the pipe. The comparison of the predicted 
pressure with the experimental data is shown in Figure 
6, which indicates that the predicted pressure is smaller 
than the measured data. 

The friction in the pipe and heat transfer between the 
gas and solid walls were neglected in the SFLOW 
results shown in Figure 6. Figure 7 shows the results 
with both friction and heat transfer effects being taken 
into account, which indicates that a better agreement 
with experimental data. This volume venting case has 
also been simulated using ISENTANK*®, which 
assumes the flow in the path is isentropic. The SFLOW 
prediction in Figure 6 is very similar to the result from 
ISENTANK with a discharge coefficient of 1.0 while 


Volume Venting Test 
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Figure 6. Comparison of the Predicted 
Pressure From SFLOW With the 
Experimental Data. (The Friction and Heat 
Transfer Were Neglected in the Prediction) 
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that in Figure 7 is similar to that from ISENTANK with 
a discharge coefficient of 0.6 


Volume Venting Test 



Figure 7. Comparison of the Predicted 
Pressure From SFLOW With the 
Experimental Data. (The Friction and Heat 
Transfer Were Accounted for in the 
Pred iction) 


RSRM NOZZLE JOINT 4 TEST 

SFLOW was also used to simulate the RSRM nozzle 
joint 4 hot fire test‘d, which was also modeled by 
Clavion''^. Figure 8 shows the location of Joint 4 in the 
RSRM nozzle and the hot fiow test fixture. The insert is 


Test Article 



Full Scale Hardware 


Figure 8. The RSRM Nozzle Joint 4 and 
Hot Flow Test Fixture 


the enlarged view of the joint 4. In this section, only 
comparisons of the SFLOW predictions with the 
measured data for configuration No. 8 at the 90 degree 
side will be discussed. The SFLOW results for other 
configurations, more results for this configuration, as 
well as the details of the model will be published later. 

A schematic of the test fixture for the RSRM nozzle 
Joint 4 hot fiow test is shown in Figure 9. There are two 
fill-tubes and bottles connected to the primary O-ring 
gland, chamfer region, as well as the secondary O-ring 
gland. The SFLOW model consists of four volumes 
(the combustion chamber, primary O-ring gland, 
chamfer region, and secondary O-ring gland) and nine 
paths (one from the combustion chamber to the primary 
O-ring gland, one from the primary O-ring gland to the 
chamfer region, one from the chamfer region to the 
secondary O-ring gland, two for the fill tubes and 
bottles connected to the primary O-ring, two for the fill 
tubes and bottles connected to the chamfer, two for the 
fill tubes and bottles connected to the secondary O- 
ring). 


Internal Test 
Fixture 
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Figure 9. Schematic of the Test Fixture for 
the RSRM Nozzle Joint 4 Hot Flow Test 

Figure 10 shows the solid grid used for the RSRM 
nozzle Joint 4 test with different colors for different 
materials. The pressure and temperature in the 
combustion chamber is specified as input for the 
SFLOW prediction (see Figure 16 for the pressure, the 
temperature is 530°R initially and increases to 4900°R 
at 0.07 sec and then remains at that high temperature). 
Impingement heat transfer is applied for several 
surfaces on the 0-rings, the inlet path, and the chamfer 
region while minor losses are applied in the paths 
where the fiow turns. Steady-state gas-only CFD 
simulations such as those discussed in Laubacher et ai 
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are used to determine which surfaces should have jci 
impingement heat transfer. 

Chamber 
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Figure 1 0. The SFLOW Grid of the Solid for 
RSRM Joint 4 Hot Test Simulation 

Isocontt)urs of the predicted temperatures of the solid at 
4.0 sec are shown in Figure I I. Initially the temperature 
is at 530°R. At 4.0 sec, the solid cells near the 
impinging surfaces as well as those near the flow path 
from the chamber to the primary 0-ring are hot while 
the temperatures at the solid cells farther away are still 
very low. 

Solid Temperatures in the Whole 

Computational Domain so.- 



Figure 11, Isocontours of SFLOW Predicted 
Solid Temperatures at 4.0 sec 


Is^jcomoLiiN of the predicted solid temperatures at l.O 
Ncj. 2.0 see. 3.0 sec. and 4 n sCw are shown in Figure 
;2. lagure ! 3. Figure 14. and lagurc 15, respectively, 
!-r part ^.t the d^Miiain inckidmg the 0-rings and 
chamicr. As time increases, more solid heats up while 
the temperatures of the solid near the impingement 
>urfaecs decrease due to heat conduction. Since the 
thermal conductivity is laigcr lor steel than carbon 
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Figure 12. Isocontours of SFLOW Predicted 
Solid at 1.0 sec 

Solid Temperatures in Part of the 
Computational Domain 
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Figure 1 3. Isocontours of SFLOW Predicted 
Solid Temperatures at 2.0 sec 
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Solid Temperatures in Part of the 
Computational Domain 



Figure 14. Isocontours of SFLOW Predicted 
Solid Temperatures at 3.0 sec 

Solid Temperatures in Part of the 
Computational Domain 



Figure 15. Isocontours of SFLOW Predicted 
Solid Temperatures at 4.0 sec 
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Fill Bottle Pressure Off 
Primary 0-ring Groove: RSRM Joint 4 



Time (sec) 

Figure 1 6. Comparison of the Predicted 
Gas Pressures With the Measured Data 


Fill Bottle Pressure Off 
Chamfer Region: RSRM Joint 4 
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Figure 1 7. Comparison of the Predicted 
Gas Pressures With the Measured Data 


Fill Bottle Pressure Off Secondary 
O-ring Groove: RSRM Joint 4 



Figure 1 8. Comparison of the Predicted Gas 
Pressures With the Measured Data 
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combustion chamber, which is an input to the SFLOW 
code. The pressures at all the fill bottles are very 
similar and the predicted values agree very well with 
the measured data. 

The predicted temperatures of the fill bottles in the 
RSRM nozzle joint 4 test are compared with the 
measured data in Figure 19. The agreement is good 
considering that the chamber temperature is about 
5,000°R while that in the fill bottles is less than 600°R. 


Fill Bottle Temperatures: RSRM Joint 4 



Time (sec) 

Figure 1 9. Comparison of the Predicted Gas 
Temperatures With the Measured Data 


Figure 20 and Figure 2 1 show the predicted solid 
temperatures at two locations, one near the primary O- 
ring (TI4 in Figure 22) and the other Just before the 
secondary O-ring (T15 in Figure 22), together with the 
measured data at the same locations. As discussed by 


Temperature Versus Time: RSRM Joint 4 



Figure 20. Comparison of the Predicted 
Temperature of the Solid Near the Primary O- 
ring (T14 in Figure 22) With the Measured Data 


Temperature Versus Time: 
RSRM Joint 4 



Figure 21. Comparison of the Predicted 
Temperature of the Solid Near the 
Secondary O-ring (T15 in Figure 22) With 
the Measured Data 



Figure 22. The Location of the Thermocouples 
In the RSRM Nozzle Joint 4 Test 

Clayton‘'\ the measured temperatures are not very 
accurate due to the large gradients, tiny gaps and brief 
time scales. On the other hand, better agreement might 
be obtained by using a full three-dimensional solid grid. 
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Total Erosion (in.) 


The predicted erosion of the secondary 0-ring is shown 
in Figure 23. The predicted total erosion after 2 sec is 
about 0.00813 in, which is very close to the measured 
value of 0.008 in. This agreement is excellent 
considering analytical modeling complexity and 
assumptions required, as well as the variability in 
measured data, measurement error due to very small 
dimensions and short time frame. . 



Figure 23. Secondary 0-ring 
Erosion Predicted by 
SFLOW In RSRM Nozzle Joint 4 Test 


As discussed above, most CFD applications in the 
RSRM joints do not include the solid region in the 
computational domain and the solid wall is assumed to 
be either adiabatic or isothermal. Figure 24 and Figure 
25 compare the SFLOW predictions of the pressures 
and gas temperatures in the fill bottles with and without 
heat transfer between the gas and solid. It indicates that, 
without heat transfer, the fill time is reduced by a factor 
of about three and the gas temperature increases to 
more than 1,400‘^R from around 580°R with heat 
transfer. For these conditions, heat transfer to the solid 
wall is a significant driver for the problem. 


Fill Bottle Pressures: RSRM Joint 4 
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Figure 24. Comparison of the 
SFLOW Predicted Gas Pressure With and 
Without Heat Transfer 
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SUMMARY AND CONCLUSIONS 

A new thermal-flow simulation code, called SFLOW, 
has been developed to model the gas dynamics, heat 
transfer, as well as 0-ring and flow path erosion inside 
the space shuttle solid rocket motor joints by 
combining SINDA/G®, a commercial thermal analyzer, 
and SHARP®, a general-purpose CFD code developed 
at Thiokol Propulsion. SHARP® was modified so that 
friction, heat transfer, mass addition, as well as minor 
losses in one-dimensional flow can be taken into 
account. The pressure, temperature and velocity of the 
combustion gas in the leak paths are calculated in 
SHARP® by solving the time-dependent Navier-Stokes 
equations while the heat conduction in the solid is 
modeled by SINDA/G®. The two codes are coupled by 
the heat flux at the solid-gas interface. 

A few test cases are presented and the results from 
SFLOW agree very well with the exact solutions or 
experimental data. These cases include Fanno flow 
where friction is important, Rayleigh flow where heat 
transfer between gas and solid is important, flow with 
mass addition due to the erosion of the solid wall, a 
transient volume venting process, as well as some 
transient one-dimensional flows with analytical 
solutions. In addition, SFLOW is applied to model the 
RSRM nozzle Joint 4 subscale hot-flow tests and the 
predicted pressures, temperatures (both gas and solid), 
and 0-ring erosions agree well with the experimental 
data. It was also found that the heat transfer between 
gas and solid has a major effect on the pressures and 
temperatures of the fill bottles in the RSRM nozzle 
joint 4 configuration No. 8 test. 
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